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The 3D one-component spin system in an external magnetic field is studied 
using the collective variables method. The integration of the partition function of 
the system over the phase space layers is performed in the approximation of the 
sextic measure density including the even and the odd powers of the variable (the 
asymmetric p 6 model). The general recurrence relations between the coefficients of 
the effective measure densities are obtained. The new functions appearing in these 
recurrence relations are given in the form of a convergent series. 



1 Introduction 

The object of investigation is the Ising model on a simple cubic lattice with 
an exponentially decreasing interaction potential. The Ising model, which 
is simple and convenient for mathematical analysis, is widely used in the 
theory of phase transitions for an analysis of properties of various magnetic 
and nonmagnetic systems (binary mixtures, lattice model of liquids, ferro- 
magnets, etc.). 

The three-dimensional (3D) Ising system in an external magnetic field 
will be studied here based on the collective variables (CV) method This 
original method can be extended to the liquid-gas system 0, multicompo- 
nent fluids ||. In this case, the method makes it possible to use a conver- 
gent series by applying the solution of the 3D Ising model in contrast to the 
asymptotic e-expansions (see, for example, @,|5)). 

The CV method allows one to calculate the partition function of the 
Ising model and to obtain not only the universal quantities but also analytic 
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expressions for nonuniversal characteristics as functions of the microscopic 
parameters of the system. The statistical description of the spin system 
behaviour is performed in real 3D space on the microscopic level. The main 
results obtained in the CV method for the 3D systems in the absence of an 
external field are presented in 

An important condition for describing the system properties by the CV 
method is the use of non-Gaussian measure densities Q. Such a measure 
density at a zero external field is represented as an exponential function of 
the CV, the argument of which contains, along with the quadratic term, 
higher even powers of the variable with the corresponding coupling con- 
stants. The confinement to the 6th power of the CV in the expression for 
a measure density corresponds to the p e model. This model provides a 
more adequate quantitative description of the 3D Ising system behaviour in 
comparison with the p 4 model (see 

In the case of a nonzero external field, the non-Gaussian measure den- 
sity includes the terms proportional to odd powers of the CV in addition 
to the terms proportional to even powers. The partition function of the 
3D one-component system in an external field is presented below in the ap- 
proximation of the sextic measure density, which involves the even and odd 
powers up to the 6th power of the variable (the asymmetric p 6 model). A 
partition function functional is similar to the Ginzburg-Landau functional 
(see, for example, flQ|] ). We calculate the partition function by successively 
integrating its expression over the short-wave or the so-called unimportant 
variables. The corresponding renormalization group (RG) transformation 
can be related to the Wilson type [11|. The general recurrence relations 
(RR) corresponding to the asymmetric p 6 model and the new special func- 
tions entering the RR are considered. 

In our earlier works (see, for example, |7],||,||]), an infinitely weak external 
field was introduced in the course of calculation of the contribution from 
the long-wave modes of the spin moment density oscillations to the 3D 
Ising system thermodynamic characteristics. In this paper, we introduce an 
external field in the Hamiltonian from the outset. Such an approach leads to 
the appearance of odd powers of the CV in the expression for the partition 
function and makes it possible to describe a lot of quantities (cumulants, 
initial coefficients of the partition function, etc.) as functions of an external 
field. 
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2 Recurrence relations for the asymmetric p 6 mo- 
del 

We consider the Ising system on a simple cubic lattice with period c. The 
Hamiltonian of the system has the form 

H = -\ E^yVi^j-^E^i- 



Here h is an external field, r;j is the distance between particles at the sites 
i and j, a\ is the operator of the z-component of spin at the ith site, having 
two eigenvalues +1 and —1. The interaction potential is an exponentially 
decreasing function 

<&(r y )=Aexp(-^), (2.2) 

where A is a constant, 6 is the radius of effective interaction. For the Fourier 
transform of the interaction potential, we use the following approximation: 

^(0)(1 -26 2 fc 2 ) at k < B', 
nk) ~ \ at B> < k < B. [2 - 6) 

Here B = tt/c is the boundary of the Brillouin half-zone, B' = is 
determined from the condition for the application of the parabolic approx- 
imation $(0)(1 - 2b 2 k 2 ) = 0, $(0) = 8irA(b/cf. At $(0) = 2dJ, 6 = 6/ = 
cj(2\fd) (J is the constant of the interaction between the nearest neigh- 
bours, d = 3 is the space dimension) for small values of the wave vectors k, 
the parabolic approximation of the Fourier transform of the exponentially 
decreasing interaction potential corresponds to the analogous approximation 
of the Fourier transform for the interaction potential between the nearest 



neighbours [12 



In the CV representation for the partition function of the 3D Ising model 
in an external field, we have 



Z = J exp 



- ]T /3^(% k P-k 
z k 



J(p) (dp) N , (2.4) 



where the summation over the wave vectors k is carried out within the first 
Brillouin zone, (3 = l/(kBT) is the inverse thermodynamic temperature, ks 
is the Boltzmann constant. The CV are introduced using the relations 
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of the type of an analytic functional for operators of spin density oscillation 
modes p^ = (x/iV) -1 J2i ci exp(— zkl) The expression for J(p) can be 
written as 



J(P) 



(2 cosh h') N I exp 



2m J>kPk + ^(-2^)"^ (n ' 2)/2 x 



■Mn 



n! 



E w ki"--^k w ^k 1 +-+k n 



ki,...,k„ 



(2.5) 



Here <5] £lH ^ is the Kronecker symbol, the variables are conjugate to 

p k . The cumulants .M n are functions of a generalized field b! = (3h: 

Mi = tanh h! = cc, 
Mi = 1 - tanh 2 /i' = y, 
M 3 = -2xy, 
Ma = -2y 2 + Ax 2 y, 
M 5 = 16xy 2 - 8x 3 y, 

M 6 = 16y 3 - 88x 2 y 2 + 16x 4 y. (2.6) 

It should be noted that (|2.5| ) is the Jacobian of the transition from the set 
of the N spin variables <j\ to the set of the CV p k at h = 0. 

We shall proceed from the expression for the partition function in the 
approximation of the asymmetric p 6 model. Putting n = 6 in ( |2.5[ ) and 
carrying out integration in fl2.4|) with respect to the variables /3 k and u; k 
with indices B' < k < B, followed by the integration with respect to the N 1 
variables w k , we obtain 



Z = (2co S hh') N 2 ( - N '~ 1)/2 e a 'o N ' |exp 



i 



fc<B' 



><PkP-k- E 



£ ^ • • • Pk ; x 



X'W-.+k, 



(dp)", 

B/B' = irV2b/c, and 
d'(k) = a' 2 -P$(k). 



where N' = Ns d , 



o ) s o 



(2.7) 



(2. 
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The coefficients 



i / 4 1 .A 1/4 
a' = lnQ(M), Q(M) = tT 1 (^o) W), 

a ' fc = te S V CfecCV), ^ = 2,3,4,5,6 (2.9) 



are functions of a field h! and a ratio b/c. In this expressions, Mi = —Mi. 
The set 7/ = {7/1, 77 2 , 7/ 3 , 7/4, 7/5, 7/g} is the set of initial arguments 

Vi - Mi {mJ s ° ' %_ 2 ^4 J ' 

% ~ 3! S ° ' 174 - 1 ' 

/ _ Ms /_4j_\ 5/4 ,_Ma / J!_\ 3 / 2 d/2 

% " 5! U 4 J ' % " 6! U 4 J S ° • (2 ' 10) 



The special functions Ci sc {r]') in (|2.9[) are combinations of the functions 
*(2/+i) sc (»/) = I(2i+i)s{v')/Ioc(v') and F 2icc (rf) = hic(v') I 'hc(rf): 

C lsc (v') = F lsc { V '), 

C 2sc ( V ') = F 2cc (v')+Ft c ( V '), 

C 3sc (7?') = F 3sc ( V ')-3F 2cc (v')Fisc(r ] ')-2Fl c (7 1 '), 

CaM) = -FaM) - 4F 3sc (r ? , )i ? is C (V) + 3F 2 2 CC (V) + 
+12^(7/0^(7/') + 6F 4 c (r/), 

Cs.cCV) = -i ? 5s C (r? , )+5i ? 4 CC (r/ , )i ? is C (r/ , ) + 10i ? 3 SC (r/')^2 CC (r? / ) + 
+20F 3sc (r/ / )Fi 2 sc (r/ / ) - 30F 2 2 cc (r/)F lsc (r/) - 

-60^(7/0^(7/0 - 24^(7/), 

C 6sc (r/0 = ^6cc(r?0 + 6^(7/0^(7/) - 15F 4cc (r/0F 2cc (7/0 - 
-30^(7/0^0(7/0 + 1OF 3 2 sc (t/0 - 

-120^(7/0^(7/0^(7/0 " 120^(7/0^(7/0 + 

+30F| cc (t/0 + 270^(7/0^(7/0 + 360F 2cc (t/)F 4 c (t/) + 
+120^(7/0- (2.11) 
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Here 



and 



oo oo 

I(2l+i)s(v') = J t 2l+1 f s (t) dt, I 2lc (r,') = J t 2l f c (t) dt, 


1 = 0,1,2,3, (2.12) 

f s (t) = sin(^i - V ' 3 t 3 + V ' 5 t 5 ) exp(-r^ 2 - r^t 4 - V ' 6 t e ), 

f c (t) = cosfai* - rj' 3 t 3 + r^ 5 ) exp(-r^ 2 ~ ^ - r^ 6 ). (2.13) 



We use the method of layer-by-layer integration of (2.7) with respect to 
the variables p^ proposed by I.R. Yukhnovskii Q. The integration begins 
from the variables p-\^ with large values of k (of the order of the Brillouin 
half-zone boundary). The phase space of the CV p^ is divided into layers 
with the division parameter s. In each rath layer (corresponding to the re- 
gion of wave vectors B n+ \ < k < B n , B n+ \ = B n /s, s > 1), the Fourier 
transform of the potential is replaced by its average value $(B n+ i, B n ). 
As a result of step-by-step calculation of the partition function, the number 
of integration variables in the expression for this quantity decreases gradu- 
ally. The partition function is then represented as a product of the partial 
partition functions of individual layers and the integral of the "smoothed" 
effective measure density. After the integration over the n + 1 layers of the 
CV space, we obtain 

Z = (2cOshti) N 2 {Nn+1 - 1)/2 QoQl---Qn[Q(Pn)] Nn+1 X 

Wt +1 \p) (dp) N ^\ (2.14) 



where N n+1 = N' S -^ n+1 \ and 



Qo = [e a '°Q(d)] N ' , Qi = [Q(P)Q(d 1 )] Nl 
Q n = [Q(P n _i)Q(d n )] Ml , 

(n)\V 4 



Q(d n ) = (4!/4 n) ) Io(h^), 
Q(P n ) = vr- 1 (s d a { fi/C A (h^jf lA hci^). (2.15) 



The components 



h (n) _ M ( 4! V /4 (n) _ d n (B n+1 ,B n ) ( 4! \ 1/2 

1 1 >) ' 2 ~ 2 >) 
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4"' = ^(4t) 3/1 . *f = i. 



3 1 




3! 1 




,(«) 
' 5 1 




5! 1 





5/4 (n) / \3/2 

■ ^-^-(Jt) < 2 ' 16 > 

of the set of basic arguments in the nth layer 

/i( n ) = {hi , /i 2 , ^3 > > ^5 > ^6 } are determined by the average value 
of the coefficient d n (k), i.e., by d n (B n+ i, B n ) = a 2 n ^ — (3&(B n +i, B n ) as well 
as the quantities at]™ , a 3 , 04 , and a 6 . The components of the set of 
intermediate arguments r/ n ) = {rj^ , 77 2 n ^ , rj^ , 774 , 775 , 77^} are functions 
of the quantities fl2.16|) . They are given by the formulas 



r,W = -(4!s M ) 1 / 4 C , 1 (/ i W)C 4 - 1/4 (/ l W), 
rj^ = (6 S d ) 1 / 2 C 2 (/ l ^)C 4 - 1/2 (/ i W), 



»7? ) - 1, 



( B ) = ( 4!g " d ) 1/4 C5(/t W )C -5/4 ( ^(n) )i 



% 5 



V M = (68 /_ )1/2 C7 6 (^)C 4 - 3/2 (fcW). (2.17) 

The special functions C[(h^) can be expressed in terms of F[(h^) = 
/K/* (n) )//o(/i (n) ): 

Ci(/iW) = Fi(fcW), 

C 2 (fr(«)) = F 2 (fcW)-if(fc< n) ), 

C 3 (/i^) = F 3 (^ n )) - 3F 2 (/i (n) )Fi(/i( n )) + 2F?(h^), 
C 4 (/t (n) ) = -F 4 (/i (n) ) + 4F 3 (/i( n ))Fi(/i( n )) + 3F$(h^) - 

-^(/i^Ff (frW) + 6if (/i^), 
C 5 (/i( n )) = -F 5 (/i( n )) + 5F4(/i (n) )Fi(/i( n )) + 10F 3 (/i( n ))F 2 (/i( n ))- 

-20F 3 (/ l (n) )F 1 2 (/ l (n) ) - 30F|(/t (n) )Fi(/iW) + 

+60F 2 (h^)F?(h^) - 24F 1 5 (/ i ( n )), 
C 6 (/t^) = F 6 (fcW) - 6F 5 (/i^)Fi(/t( n )) - 15F4(/i (n) )F 2 (^ n )) + 

+30F 4 (/i (n) )F 2 (/iW) - 10F 2 (/i( n )) + 
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Here 



+l20F 3 (h^)F 2 (h ( - n) )F 1 ( y h^) - 120F 3 (/i( n ))F 1 3 (/iM) + 
+30Ff(/i( n )) - 270F|(/ l ( n ))F 1 2 (/ l ( n )) + 
+360F 2 (/ i ( n ))F 1 4 (/ i ( n )) - 120Ff(hF>). (2.18) 



Ii(hW)= Jt l f^(t)dt, 1 = 0,1,2,3,4,5,6, 



(2.19) 



and 



f^(t) = eM-hT't ~ hT'f ~ hf ! t A - hT't 11 - hf't* - h^>e). (2.20) 

The effective sextic measure density of the (n + l)th block structure 
Wg" - " 1 " 1 -* (p) is written as 



»+2 



(n),3 lW,4 j»+5 



(n+1) 



(p) = exp 



fc<_B n+ i 



6 „(n+l) 



Z=3 ^ n +l fci,...,fci<B„+i 



1/2-1 



Pki ■ ■ ■ Pk ; X 



X*k 1+ ...+k, 



(2.21) 



a ; are the renor- 



where B n+1 = B's~^ n+1 \ d n+1 (k) = a { 2 n+1) - 0$(k), 
malized values of the coefficients a[ after the integration over the n + 1 layers 
of the phase space of the CV. 

The coefficients of the sextic measure densities of the (n + l)th and nth 
block structures are connected through the following general RR: 



-^- 2 + (4 n) / /4 ^ (n) )' 



I = 1,2,3,4,5,6. (2.22) 



We have introduced new designations here: 



u 



in) _ 



s a 



(n) 
1 > 



4 n) +q = s 2n d n (B n+1 ,B n ), 



4 n) = s 2n d n (0), 



u 



(n) _ mn (n) 



m = 3, 4, 5, 6. 



(2.23) 



The quantity q = qy3$(0) determines the average value of the Fourier trans- 
form of the potential in the nth layer (3§(B n+ i,B n ) = /34>(0) - q/s 2n (q 
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corresponds to the average value of k 2 on the interval (1/s, 1]). The func- 
tions appearing in the RR ( gjp 

can be defined by the expressions 
Yl (h^) = s- d / A C ls M n) )Cl l, \h^), 

k = 2,3,4,5,6. (2.24) 

The functions C; sc (r/^) have the form similar to C; sc (r/) (|2.1l| ) under con- 
dition that the set r/ and the functions f s (t), f c (t) ( |2.13| ) should be replaced 
by i/ n ' and 

F,(t) = sm(r,[ n) t + r/^t 3 + V ^t 5 ) expi-^t 2 - V ^t 4 - r,^t% 

Fc (t) = cosfa^t + ^k 3 + ^ n) i 5 ) exp(-^ n) t 2 - i^V - ^ n) t 6 ), (2.25) 



respectively. It should be noted that the formula for I^djnf) ( see ( 2.12 )) 
allows us to express the function Iodv^) entering ( |2.15| ) in the same way. 
The function C 4 (h^) is given in ( p7i8|) . 

The obtained expressions make it possible to investigate the properties 
of the 3D one-component spin system in the vicinity of the critical point 
(T = T c , h = 0). In the absence of an external field (h = 0), we have 
= a^/ = a^ 1 = and relations ( 2.22 ) are indentical with the RR 
presented in [13|. The behaviour of the system can also be studied at T = T c 
and h — ► 0. 

The components of the set of intermediate arguments r/ n ) and the func- 
tions Ci sc (rj^), Ci(h^), Yi(h^) are approximated by power series in de- 

(n) 

viations of the basic arguments h\ from their values at a fixed point h i . 
Taking into account the linear deviations, we find the following forms for 
these series: 

A x = A m [l + A a {h\-ht ) )+A l2 {hl-h^ ) ) + A m (hl-hf ) )+ 

+ A lb {h%-hf ) ) + A l s(h%-hf } )\, 1 = 1,2,3,4,5,6. (2.26) 

Here A\ are the above-mentioned intermediate arguments or functions. The 
coefficients An are ultimately functions of the basic arguments at a fixed 
point. It should be emphasized that we necessarily arrive at asymptotic 
expansions when the Gaussian density of measure is taken as the basis den- 
sity. The use of the sextic non-Gaussian measure density permits to ob- 
tain the rapidly converging series (HH). Having (B) , we can write and 
study the approximate RR. On the basis of ( p.l4j) and RR, the free energy 
F = — ksTlnZ can be calculated. 
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3 Discussion and conclusions 



The results obtained for 3D Ising model in an external field can be applied 
to the description of liquid-gas critical points of both a one-component fluid 
P,|l4| and a binary fluid mixture (see, for example, fl5p . The CV method 
with a reference system (RS) is used for these systems. The expression 
for the partition function contains the even and odd powers of the variable 
with the corresponding cumulants and is similar to ( |2.4j ), (|2.5|). Detailed 
investigation of the properties of RS cumulants in Q makes it possible to 
transform the grand partition function into a functional form defined on the 
effective block lattice. The functional corresponds to the partition function 
of the Ising model in an external field. A new point in problem of the 
liquid-gas critical point as compared to the case of the Ising model is the 
dependence of the critical temperature and all coefficients on the density and 
chemical potential. The latter is equivalent to the insertion of a constant 
external field into the Ising model. To describe the effects connected with the 
asymmetry for the surface of coexistence of the phases, one should compute 
the partition function using the sextic measure density [^, |l6| . 

In this paper, the 3D one-component spin system in a nonzero external 
field is investigated by the CV method in the asymmetric p e model ap- 
proximation. An initial expression for the partition function of the system 
is constructed in the form of a functional with explicitly known coefficient 
functions (see (|2.7D ). The partition function is integrated over the layers 
of the CV phase space. The main feature is the integration of short-wave 
spin density oscillation modes, which is generally done without using pertur- 
bation theory. The short-wave modes are characterized by the presence of 
the RG symmetry and are described by the non-Gaussian measure density. 
The general RR ( 2.22| ) for the coefficients of the even and the odd powers 
of the variable in the two adjacent sextic measure densities are found. The 
new special functions appearing in the construction of the theory using the 
asymmetric p 6 model are considered. These functions entering the RR are 
represented in the form of a power series in small deviations of the basic 
arguments from their values at a fixed point. Representing the RR in the 
form of a nonasymptotic series relates to rejecting the traditional use of 
perturbation theory, which is based on the Gaussian measure density. 

The expressions under investigation allow us correctly to trace the de- 
pendence of the results on a field within the framework of the CV method. 
As a consequence of the presence of an external field h in the Hamiltonian 
(|2.lD, we obtain odd powers of the variable in the expressions (|2.5D for J(p) 
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and ( p.7| ) for the partition function Z. The cumulants M n in Q2.5| ) are com- 
binations of the quantities x = tanh/i' and y = 1 — tanh 2 h', where h' = f3h 
is a generalized field. The initial coefficients a\ (2.9) of the partition function 
are determined by these cumulants and are depended by sq on the micro- 
scopic parameters of the system (the effective interaction radius b and the 
lattice constant c). The cumulants Mi, M3, .M5 as well as the coefficients 
a'i, a' 3 , a' 5 are odd functions of an external field. The quantities M2, M4, 
M.q and a' , a' 2 , a' 4 , a' 6 are even functions of a field. The RR ( |2.22 ) make 
it possible to determine the elements of the block Hamiltonians from the 
initial data. As follows from the obtained expressions, the cumulants A4i, 
M3, M-5 and coefficients a[, a' 3 , a' 5 , ... , a{ , a 3 , 05 vanish at h = 0. 
Then the RR ( [2.22 ) for the asymmetric model turn into the RR for the 
p 6 model with even powers only [13|. In the case when h 7^ 0, the measure 
density involves odd powers of the variable p-^ in addition to the even pow- 
ers. Although the CV method as well as Wilson's approach exploit the RG 
ideas, it is based on the use of non-Gaussian measure densities. This allows 
one to obtain a qualitatively new form of the RR between the coefficients of 
the block Hamiltonians. In the limiting case (corresponding to the Gaussian 
basis measure density; h = 0) these RR reduce to the Wilson RR [11|. As 
was shown in [17], while this limiting case does not allow one to perform the 
calculation of the expression for the free energy of the system, it provides 
reliable results for the critical exponents of thermodynamic characteristics. 
The solutions of the RR under consideration can be used for calculating 
the thermodynamic functions of the system and for deriving the equation of 
state using the method proposed in [111] . 
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